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whittaker limits of difference spherical functions 3 

0. Introduction 

The main aim of this paper is to introduce global g— Whittaker func- 
tions as the hmit t —>■ oi the (renormahzed) generahzed symmetric 
spherical functions constructed in [C5] for arbitrary reduced root sys- 
tems (see [Sto] in the C^C— case). This work is inspired by [GLOl] 
and [GL02], though our approach is different. For instance, we ob- 
tain a g— version of the classical Shintani-Casselman-Shalika formula 
[Shi, CS] via the g— Mehta -Macdonald integral in the Jackson setting. 
The Shintani-type formulas (in the case of GLn) play an important role 
in [GLOl, GL02], but the g— Gauss integrals are not considered there 
as well as globally-defined g— Whittaker functions. We use these for- 
mulas to obtain a g, t— generalization of the Harish- Chandra asymptotic 
formula for the classical spherical function. 

0.1. Results and applications. The key observation is that the def- 
inition of the symmetric g, t— spherical functions from [C5] is compat- 
ible with taking the Whittaker limit and results in globally-defined 
g— Whittaker functions. The definition from [C5] is based on the g— 
Mehta- Macdonald integrals calculated there for the constant term 
functional, i.e., in the setting of Laurent series. In this paper, we mainly 
treat the spherical functions as global ones, analytic or meromorphic. 

The g— Whittaker functions are solutions of the g— Toda eigenvalue 
problem and are expected to have important applications in math- 
ematics and physics, including the Langlands program. Concerning 
the latter and relations to the affine flag varieties, see, for instance, 
[GiL, BF, Ion2]. The g— Shintani-Casselman-Shalika formula gives a 
(relatively simple) example of the Langlands correspondence. The 
affine Toda lattice provides another link; it is (presumably) dual to 
the g— Toda lattice in the sense of [KL] (via the monodromy map). 

The coefficients of the expansion of our Whittaker function are essen- 
tially polynomials in terms of g with non-negative integral coefficients. 
It can be verified using the intertwining operators or via the relation to 
the Demazure characters (we do not discuss it in this paper). This fact 
is of obvious importance for the "categorization" of the g— Whittaker 
function and its geometric applications. The "g— integrality" has no 
known counterpart in the general g, t— theory (with a reservation con- 
cerning the stable GL— case); one of the parameters, g or t, has to be 
eliminated or expressed in terms of the remaining one. However, the 
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g, t— spherical functions are more convenient to deal with in many other 
aspects, including the analytic theory. 

The cornerstone of their theory is the duality based on the DAHA- 
Fourier transform; see [C8] and[C9]. It is a special feature of the gen- 
eral q, t— theory, missing for the g— Whittaker functions and g— Hermite 
polynomials (t — 0), in the Harish- Chandra theory (g 1) and in 
the p— adic limit (g oo). A specific problem with the Whittaker 
limiting procedure among other degenerations is that it destroys the 
VT— invariance; it gives another reason for treating Whittaker functions 
as limits of the spherical functions rather than for creating their intrin- 
sic theory. On the other hand, the Whittaker functions satisfy quite a 
few identities that are not present in the g, t— theory. These formulas, 
the g— integrality of the coefficients and various applications obviously 
make the g— Whittaker functions an important independent direction, 
which requires developing specific methods. 

At the end of the paper, we outline the approach to the global sphe- 
rical and Whittaker functions via the harmonic analysis. Our formulas 
for these functions are actually equivalent to certain fundamental prop- 
erties of the corresponding integral transforms in the space of Laurent 
polynomials multiplied by the Gaussian. The latter space is the sim- 
plest and the most natural choice here, but the same functions can 
serve other algebraic and analytic situations. The harmonic analysis 
direction seems very promising. For instance, the existence of the g— 
Whittaker limit of the global spherical function appears a boundary 
case of the general theory of growth estimates for the g, t— spherical 
function in terms of x, A (the spectral parameter) and t = . 

0.2. Growth estimates. Provided that 3?(x), 3?(A) are inside the pos- 
itive Weyl chamber (the walls must be avoided), the global spherical 
function for < g < 1 approaches asymptotically in the limit of large 



for the theta-series and g— Gamma function associated with a given 
root system R; t = q^, Pk = kp in the simply-laced case, CT being the 
constant term of the celebrated (symmetric) Macdonald function. No 
inequalities for k are necessary but one must avoid the values where 
the polynomial representation of DAHA becomes non-semisimple. 
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Up to a periodic function, the x— dependence of this function is 
this theorem is an exact g, t— analog of the Harish- Chandra 
formula [HC] describing the asymptotic behavior of the classical sphe- 
rical function in terms of the c— function. In the Whittaker limit, pk is 
omitted and 3?(a;) must be taken from — C+ (the Whittaker function it 
is not M^— invariant with respect to x). 

It is one of the major results of this paper, which seems a beginning 
of fruitful analytic g— theory. 

0.3. Our approach. It is different from that of [GLOl, GL02] (and 
we deal with arbitrary reduced root systems). The technique of the 
Gaussians is the key to introduce the global g— Whittaker function and 
prove the Shintani-type formulas. The g— Whittaker function is mainly 
treated in [GLOl, GL02] as a discrete function on the weight lattice 
for GLn satisfying the g— Toda system of difference equations. 

The space of all solutions is, generally, dimensional over the field 
of periodic functions, playing the role of constants in the difference 
theory; upon the restriction to the weight lattice it is dimensional 
over C. Ghoosing the "right" Whittaker function in this space re- 
quires certain growth conditions; using the VT— symmetric dependence 
on the spectral parameters gives another approach. There is no in- 
trinsic definition of the g— Whittaker function so far, but our formula 
and the growth conditions we establish clarify what can be expected. 
First, only positive powers appear in its Laurent series expansion (af- 
ter dropping the Gaussians). Second, our x— asymptotic formula for the 
g— Whittaker function inside the negative Weyl chamber is sufficient to 
fix it uniquely. 

We note that in the differential setting, the spherical and Whit- 
taker functions can be uniquely determined from the eigenvalue prob- 
lem (subject to the TV— invariance for the spherical function and certain 
growth conditions in the Whittaker case). It simplifies the starting def- 
initions. However the difference theory is more universal and, remark- 
ably, has important algebraic and analytic advantages. The self-duality 
of the DAHA- Fourier transform and the technique of the Gaussians are 
the key; these are special features of the g, t— setting and are mainly ab- 
sent in the trigonometric-differential and p— adic cases. In this respect, 
the g, t— theory is somewhat similar to the rational-differential theory 
of (multi-variable) Bessel functions. 
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0.4. Difference spherical functions. The global nonsymmetric and 
symmetric q, t— spherical functions were defined in [C5] and then in [Sto] 
(the C'^C— case) as the reproducing kernels of the Fourier transform of 
the standard polynomial representation twisted by the Gaussian. In 
this approach, the spherical function is determined uniquely (the Mac- 
donald eigenvalue problem fixes it only up to periodic factors). Using 
the Gaussians, among other things, provides the global convergence. 
These construction appeared compatible with the Whittaker limit. 

The Gaussians play the key role in our approach to the Shintani- 
Casselman-Shalika formula. In the g, t— setting, it becomes the Mehta 
-Macdonald formula in the Jackson case from [C5], where a special 
vector, —pk, is taken as the origin of the Jackson summation. 

Developing this direction, we conclude the paper with the Jackson- 
Gauss integrals for the global spherical and Whittaker functions; such 
formulas were given only for Macdonald polynomials in [C5] . These for- 
mulas seem an important step toward systematic difference harmonic 
analysis, although the case of the real integration is still beyond the ex- 
isting theory. Now, with the q, t—Harish- Chandra asymptotic formulas 
from this paper, it seems that there are no obstacles for developing the 
real integration theory generalizing the classical "non-compact" case. 

Conceptually, as it was observed in [GL02], the g— variant of the 
Shintani-Casselman-Shalika formula is nothing but the duality formula 
for the Macdonald polynomials from [C3] considered upon the limit 
t ^ 0. However, establishing exact relations is, generally, a subtle 
problem. The Shintani-type formulas play the major role in the paper, 
including the growth estimates. 

This interpretation gives evidence that the DAHA- Fourier transform 
is connected with the (local quantum) geometric Langlands correspon- 
dence. Generally, the DAHA— localization functor, which includes the 
modular transformation q i— s> q', is expected to play its role in the quan- 
tum geometric Langlands correspondence; the DAHA-Fourier trans- 
form is likely to be one of its ingredients. 

We note that DAHA leads to a theory that is a priori more general 
than the one needed for the (local) quantum Langlands correspondence 
because it contains an extra parameter t. However, there is growing ev- 
idence that the general q, t— DAHA appear in the Langlands program. 
It makes important the exact relations between the g, t— spherical func- 
tions and g— Whittaker ones (which are already a part of the Langlands 
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program). We expect this paper to trigger interesting new develop- 
ments. 

It is worth mentioning that the approach to spherical functions via 
the Fourier transform depends on the choice of the corresponding rep- 
resentation of the double afiine Hecke algebra. Technically, the choice 
of this space influences only the normalization; spherical function are 
defined up to periodic factors. However, the analytic properties of the 
g, t— spherical function, exact factors in the Shintani-type formulas and 
other similar features reflect the properties of the considered represen- 
tation (equivalently, the choice of the normalization). 

For instance, if the Gaussian is interpreted as a theta-function, then 
the corresponding spherical function is meromorphic but not analytic. 
Treating the Gaussian as (not as a Laurent series), i.e., using a 
somewhat different analytic setting, leads to the q, t— spherical functions 
analytic everywhere, but not single-valued in terms of q^. If the Gaus- 
sians are omitted in this definition, i.e., the DAHA- Fourier transform 
acts from the polynomial representation to the space of delta-functions, 
then the corresponding spherical function will become a generalized 
function. A general problem is in finding a representation that ensures 
the best analytic properties of the reproducing kernel; if ^ 1, then 
the polynomial representation times the Gaussian is the one. 

0.5. The setting of the paper. Only the symmetric theory will be 
considered in this work; the (truly) nonsymmetric g— Whittaker func- 
tion can be defined as certain limits of the nonsymmetric global spheri- 
cal function, but the construction becomes more involved and will be a 
subject of the next work(s). Nevertheless, we begin the paper with the 
account of the nonsymmetric Macdonald polynomials including their 
(straight) degeneration as t ^ 0, which is closely related to the De- 
mazure characters of irreducible affine Lie algebras; see [San, lonl]. 
We mainly need the formulas in terms of the intertwining operators 
to justify some of our claims and estimates; the intertwiners can be 
naturally defined only in the nonsymmetric theory. 

We mention that the Macdonald symmetric polynomials considered 
under the limit t — > generalize the classical g— Hermite polynomials, 
so the main result of the paper is in establishing the formula for the q— 
Whittaker function in terms of multi-variable q—Hermite polynomials. 

In the theory of nonsymmetric Whittaker functions (it is beyond 
this paper and not completed so far), the nonsymmetric g— Whittaker 
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function become a generating function for all Demazure characters^ 
not only the ones for ant i- dominant weights. To be exact, g— Hermite 
polynomials appear here instead of the Demazure characters (there is 
a direct link). However, this interpretation requires new technique of 
VT— spinors, and the analytic aspects are not clear at the moment. The 
appearance of all Demazure characters can clarify the role of Whit- 
taker functions in the Kac-Moody theory and may have connections to 
[GiL] (quantum i^'— theory of affine flag varieties) and to questions and 
conjectures from [BF] concerning the /C— theory of affine flag varieties. 

We note that there are two possible setups in the DAHA theory for 
the non-simply-laced root systems, which correspond to two possible 
choices of the affine extension. In this paper, we introduce the affine 
root system using ao = [—'^i 1] in terms of the maximal short root d 
(the so-called twisted case). The conjugation by the Gaussian and the 
Fourier transform preserve the double affine Hecke algebra in this setup. 
By the way, it is exactly the case when a relation to the Demazure 
characters can be established according to [lonl]. Theorem 1. 

The case of the "standard" non-twisted affine root system with = 
[—6, 1] for the maximal long root 6 is analogous, although the Fourier 
transform acts from the double affine Hecke algebra to its dual in the 
B, C— cases. This setting is expected to be related to the geometric 
Langlands correspondence (cf. [C9]). 

Technically, the switch to the standard DAHA can be achieved by 
changing the action of Tq in the polynomial representation. This change 
influences the relations of Tq with the X— operators (indexed by the 
weights). The F— operators become labeled by the coweights for such 
choice of Tq; they are labeled by the weights in this paper. However, 
this transformation is far from being direct at level of difference Mehta- 
Macdonald formulas we need for the theory of difference spherical and 
Whittaker functions. 

Acknowledgements. The author is thankful to D. Kazhdan for 
alerting me to the works of Gerasimov et. al and for our various con- 
versations on the Whittaker functions and the Langlands correspon- 
dence. I indebted to D. Gaitsgory for the discussion of the quantum 
geometric Langlands duality. Special thanks go to A. Gerasimov for 
his explanations of the results of [GLOl, GL02], which influenced this 
paper a great deal. I am very grateful to E. Opdam and J. Stokman 
for reading the paper and suggesting various improvements. 
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1. Double Hecke algebra 

Let R = {a} C M" be a root system of type A, B, ...,F,G with re- 
spect to a euchdean form [z, z') on M" 3 z, z' , W the Weyl group gen- 
erated by the reflections Sa, R+ the set of positive roots = —R+) 
corresponding to fixed simple roots ai,...,a„, F the Dynkin diagram 
with {oi, 1 < z < 77,} as the vertices. Accordingly, 

= {ct^ = 2a/ {a, a)}. 
The root lattice and the weight lattice are: 

Q = ®ti^ai CP = ©ILi^^i, 

where {ui} are fundamental weights: = 6ij for the simple 

coroots ttj^. Replacing Z by Z-t = {m G Z, ±m > 0} we obtain Q±, P±. 
Here and further see [B]. 

The form will be normalized by the condition (a, a) = 2 for the short 
roots in this paper. Thus, 

= (a, a)/2 can be either 1, or {1, 2}, or {1, 3}. 
This normalization leads to the inclusions Q C Q^,P C P^, where 
is defined to be generated by the fundamental coweights {uj^} dual 
to {ai}. 

We set z/j = Uoii, J^r = {^a,C( G R} and 
(1.1) p„ = (1/2) a = cUj, where a G P4-, z/ G z^r. 

Note that (p^, a{) = 1 for z/j = z/. 

1.1. AfRne Weyl group. The vectors 5 = [a, z/^j] G R" x M c M''+^ 
for a G P, j G Z form the ajfine root system R D R {z ^ are 

identified with [2^,0]). We add ao = to the simple roots for 

the maximal short root G P+. It is also the maximal positive coroot 
because of the choice of normalization. 

The corresponding set P+ of positive roots equals P+U{[a, z^aj], a G 
P, j > 0}. Indeed, any positive affine root [a, z/^j] is a linear combi- 
nations with non- negative integral coefficients of {a^, < i < n}. 

We complete the Dynkin diagram F of P by (by to be more 
exact); it is called ajfine Dynkin diagram F. One can obtain it from 
the completed Dynkin diagram from [B] for the dual system P^ by 
reversing all arrows. 



10 IVAN CHEREDNIK 

The set of the indices of the images of ao by all the automorphisms 
of f will be denoted by O (O = {0} for Es,Fi,G2). Let O' = {r e 
0,r ^ 0}. The elements Ur for r & O' are the so-called minuscule 
weights: {uJr,a'^) < 1 for a G R+. 

Given a = [a, Uaj] G -R, b ^ P, let 

(1.2) S5(i) = i'-(^,«^)5, b'{I) = [zX-{z,b)] 

hil = [z,C] _ _ 

The affine Weyl group W is generated by all (we write W = 
{sa, ot G R+))- One can take the simple reflections Sj = s^. (0 < z < n) 
as its generators and introduce the corresponding notion of the length. 
This group is the semidirect product W^Q' of its subgroups W = 
{Sa, a G -R+) and Q' = {a', a G Q}, where 

(1.3) a = SaS[a,u^]= S[^a,u^]Sa for « G -R. 

The extended Weyl group W generated by W and P' (instead of 
Q') is isomorphic to VTxP': 

(1.4) {wb'){[z,C]) = [w{z)X-{z,b)] for weW,beB. 

From now on, b and b', P and P' will be identified. 

Given b E P+, let Wq be the longest element in the subgroup Wq C W 
of the elements preserving b. This subgroup is generated by simple 
reflections. We set 

(1.5) Ub = WqWq eW, VTfe = b{ubY^ G W, Ui = u^^^-Ki = vr^,, 

where wq is the longest element in PF, 1 < i < n. 

The elements vr^ = vr^^ ,r & O' and ttq = id leave F invariant and 
form a group denoted by 11, which is isomorphic to P/Q by the natural 
projection {ur ^ ii'r}- As to {n^}, they preserve the set {—'d, ai, i > 0}. 
The relations Trr{ao) = ctr = (^r)~^(— "*?) distinguish the indices r G O'. 
Moreover, 

(1.6) W = Ili><W, where 7rrSi7r~^ = sj if 7rr{ai) = aj, 0<j<n. 

We will need the following affine action of on z G M"': 
{wb){{z)) = w{b + z), w eW,be P, 

(1.7) ssdz)) = z - {{z,a'^) +j)a, a = [a,Uaj] G R. 
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For instance, (6w)((0)) = b for any w & W. The relation to the above 
action is given in terms of the affine pairing {[z^ /], z' + d) = (2;, z')+l : 

(1.8) {w{\z,l]),wiz'))+d) = {[zj],z' + d) for w eW, 
where we treat d formally. 

1.2. The length on W. Setting w = HrW G IV, vr^ G H, u; e H^, the 
length l{w) is by definition the length of the reduced decomposition 
w = Siy.Si^Si^ in terms of the simple reflections Sj,0 < i < n. The 
number of Sj in this decomposition such that z/j = z/ is denoted by 
lu{w). 

The length can be also defined as the cardinality |A('u;)| of the \—set 
of w : 

(1.9) \[w) =R+n w'\R_) = {ae R+, w{a) E R.}, weW. 
Alternatively, 

(1.10) X{w) = UuX^{w), Xu{w) = {a e X{w),u{a) = u}. 

The coincidence with the previous definition is based on the equiva- 
lence of the length equality 

(1.11) (a) l^iwu) = liyiw) + l,y{u) for w,ueW 
and the cocycle relation 

(1.12) (6) X,{wu) = X,{u)Uu-\X,{w)), 
which, in its turn, is equivalent to the positivity condition 

(1.13) (c) u-\X,{w))cR+ 
and is also equivalent to the embedding condition 

(1.14) (d) X,{u) C X,{w). 

See, e.g., [C4, C8] and also [B, Hu]. Applying (1.12) to the reduced 
decomposition w = UrSii ■ ■ ■ Si^Si^, 

X{w) = { 5' = w'^SiXoiiX . . . , 5^ = Si^Si^{ai^), 

(1.15) 5^ = 5^,(0^2), 5^ = Oil }. 
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1.3. Reduction modulo W. It generalizes the construction of the 
elements Hb for b E P^; see [C4] or [C8]. 

Proposition 1.1. Given b G P, there exists a unique decomposition 
b = TTbUb, Ub E W satisfying one of the following equivalent conditions: 

(i) l{'Kb)+l{ub) = l{b) and l{ub) is the greatest possible, 

(ii) \{TXb)nR = 0. 

The latter condition implies that l{7Tb) + l{w) = l{'Kbw) for any 

w E W. Besides, the relation Ubip) = b-. E P- = —P+ holds, which, 
in its turn, determines Ub uniquely if one of the following equivalent 
conditions is imposed: 

(Hi) l{ub) is the smallest possible, 

(iv) if a E \{ub) then {a,b) ^ 0. 

□ 

Condition (ii) readily gives a complete description of the set vrp = 
{tt;,, E P}, namely, only [a < 0, v^j > 0] can appear in X^iib). 

Explicitly, 

(1.16) \{b) = {5 > 0, (6, a"^) > J > if aE R+, 

{b.a"^) > J > if aE R.}, 

(1.17) Xi-Kb) = {5 > 0, a e (6„, a"^) > j >0 if u^\a) E R+, 

{b-, a"^) > J > if u~\a) E 

For instance, l{b) = l{b^) = -2(p^,6_) for 2p^ = Ea>o"^- 

The element b- = Ub{b) is a unique element from P_ that belongs 
to the orbit W{b). Thus the equality c_ = 6_ means that 6, c belong 

to the same orbit. We will also use 6+ = Wo(fe-), a unique element in 
W{b) n P+. In terms of Tr^, 

Ub-Kb = b_, TTbUb = 6+. 

Note that l{TTbw) = l^iTb) + l{w) for all 6 G P, w E W. For instance, 

(1.18) l{b.w) = /(6„) + l{w), l{wb+) = l(b+) + l{w), 
l{ub'Kbw) = l{ub) + l{Tib) + l{w) for b E P, w E W. 



Partial ordering on P. It is necessary in the theory of nonsym- 
metric polynomials. See [Op, M3]. This ordering was also used in 
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[C2] in the process of calculating the coefficients of operators. The 
definition is as follows: 

(1.19) b<c,c>b for b,ceP if c-beQ+, 

(1.20) b^c^chb if 6_ < c_ or {6_ = c_ and b < c}. 

Recall that 6_ = c_ means that b, c belong to the same IV— orbit. We 
write <,>,-<, :^ respectively if 6 7^ c. 
The following sets 

(r{b) = {ceP,cy b}, a,{b) = {c e P, c ^ 6}, 

(1.21) a_(6)=a(6„), a+(6) = (t,(6+) = {c G P, c„ > 

are convex. By convex, we mean that if c, d = c + ra G cr for a G 
P+, r G Z+, then 

(1.22) {c, c + a, c + (r - l)a, d} C a. 



1.4. More notations. By m, we denote the least natural number such 

that (P, P) = (l/m)Z. Thus m = 2 for -D2A:, m = 1 for P2fc and Cfc, 

otherwise m = |n|. 

We will need to include the case t = in our definition, which 

requires minor deviations from the definitions of [C8],[C4] and other 

1/2 

author's papers. Namely, we multiply all Tj there by and change 
the formulas correspondingly. 

The double affine Hecke algebra depends on the parameters g, tj^, v & 
{fa}- It will be defined over the ring Q[g'''^/™', t^] formed by polynomials 
in terms of q^'^/^ and {ty]. We will also use a greater ring 

Q'^^^ ^ {c G Q(g^^/",t^) I c is well defined when U = 0}, 

which is a subring of the field of fractions of Q[g^^/"\ t^] 
We set 

(1.23) where a = [a, Uaj] G P, < i < n. 

It will be convenient to use the parameters {k^} together with {ti,}, 
setting 

ta = tu = qt for ^ = ^a, and pk = (1/2) y^fcaQ. 

Q>0 
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Note that {pk.a^) = h = = {{PkY.Oii) for i > 0; (p^)^ = 
YIi^Ap^Y ■ Using that Wo{pk) = —pk, "we obtain that {pk, —'Wo{b)) = 

{pk,b). For instance, {pk,b+) = —{pk,b-), where 6+ = wo{b^) (see 
above). 

By we mean H^gva ^^^^ a e R, {p^Y = p^/u, and 

this product contains only integral powers of tght and ting. 
For pairwise commutative Xi, . . . , Xn, 

n 

(1.24) X~, = Hx^q^ if b= [b,j], w{X~) = X^^^y 

i=l 

1 — 

where 6 = koJi e P, j e — Z, w eW. 
^-^ m 

i=l 

For instance, Xq = Xa^ = qX^^. 

We set (6, c) = {b, c) ignoring the affine extensions in this pairing. 

1.5. Main definition. We note that vr""^ is vr^* and is Ur* for 
r* E O , Ur = rc^^Ur- The reflection * is induced by an involution of 
the nonaffine Dynkin diagram F. 

Definition 1.2. The double affine Hecke algebra 7H is generated over 
Q[g±i/™^t^j iiy tfiQ elements {Tj, < ? < n], pairwise commutative 
{Xft, b G P} satisfying (1.24), o.i^'d the group 11, where the following 
relations are imposed: 

(0) {Ti-ti){Ti + l) = 0, < z < ri; 

(1) TiTjTi... = TjTiTj..., rriij factors on each side; 
(a) TXrTi'K''^ = Tj if 7rr{ai) = aj; 

(i%%)T,X^ = X^X-l{tiT-^} tf (6,ar) = l, < z < n; 

(iv) TiX, = X,Ti if {b, a^) = for 0<t<n; 

(v) vr.X.vr-i = = X^-.^,^q^^^*^'\ r e O' . 

Here and further the brackets {•} will be used to show explicitly 
the elements from t— localization of 7H that belong to 7H , i.e., those 
that do not involve t~^ and other negative powers of ty. It is not a 
new definition, but can help the readers to see which operators are 
actually from 7H \, quite a few (transitional) operators will involve t~^. 
We will postpone with the independent theory of nil-DAHA, the limit 
of 7H as t — * 0, till the next paper (s). In this paper, we use the 
standard theory of DAHA when convenient (which requires t~^). The 
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key examples are the elements {ty^Tf^} which do belong to Tti thanks 
to the renormalization. 

One can rewrite (iii,iv) as in [L]): 

(1.25) T,X, - X,^(,)T, = - l) ^y^lf' , < z < n. 

Given w e W, r E O, the product 

/ I 

(1.26) Tt,,{s = 7irY\_Tik, where w = Y\_Si^J = K^)^ 

k=l k=l 

does not depend on the choice of the reduced decomposition (because 
Ti satisfy the same "braid" relations as Si do). Moreover, 

(1.27) T^Tyj = TisQ wheneYei l{vw) = l{v) + l{w) for v.w^W. 
In particular, we arrive at the pairwise commutative elements: 

n n 

(1.28) n = q^^+-^^P^^ J] if = 5Z UuJ^ G P, ^ T^^,b e P. 

1=1 1=1 

The factors here are needed to make them from ftC ; 6+ is a unique 
element in W{b) n P+. Note that YbY-b = g^C^+'pfe). 

Generally, if we replace Si by Tj or Tf^ in any reduced decomposition 
of w; G W, then such product belongs to 71rC upon the multiplication 
by the product of ti corresponding to the terms T~^. 

The relations dual to (iii,iv) hold (for i > only): 

{t,T-i}n = n,(fe)T, if {b,a^) = l, 

(1.29) T^n = if {b,a'^) = 0, 1 < i < n. 

The counterpart of (1.25) is as follows: 

(1.30) r,n - n,(6)T. = iU - 1) , ' 1 ^ ^ ^ ^' 

where 6' = 6,-& respectively for long, short (it is the only root in the 
intersection W{ai) fl P+). 

Here and below we use that given b G P, replacing all T^^^ by 
tf^ in the product of (1.28) for Yj, results in the t— power g^^/'*'^) = 
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In the standard DAHA theory, q '^^+^Pk)Y^ for any h can be repre- 
sented as the product Tir{tJ^^'^T^'^) ■ ■ ■ {tf^^'^T^^) for a given reduced 
decomposition h = i^rSii ■ ■ ■ and proper choice of {±}. The number 
of terms is / = = 2(p^, Only positive powers T^^ will appear 
in this product when b G P+. The total number of the terms T^^ with 
Ui = u in this product equals 2((pj,)^, 

2. Polynomial representation 

From now on, we will switch from Tti to its intermediate subalgebra 
TH!' C ?H with P replaced by a lattice B between Q and P (see [C7]). 
Accordingly, 11 is changed to the preimage of B/Q in n. Generally, 
there can be two different lattices Bx and By for X and Y. We consider 
only Bx = B = By in the paper; respectively, a,b E B in Xa, Yf,. 

We also set W'' = B ■ W C W, and replace m by the least m E N 
such that m{B^ B) C Z in the definition of the Q'^f 

Note that TH}' and the polynomial representations (and their rational 
and trigonometric degenerations) are actually defined over Z extended 
by the parameters of DAHA. However the ring Q'^ ^ will be sufficient in 
this paper. 

The Demazure-Lusztig operators follows: 

(2.1) Ti = t,Si + {ti - 1)(X„^ - \Y\si - 1), < « < n; 

they obviously preserve Q[g, We note that only the formula for 

To involves g: 

= toso + (to - l)(Xo - l)"^(so - 1), where 

(2.2) Xo = qX-\ So{X,) = X.X;^"'^^ q^^'^l «o = h^, !]• 

The map sending Tj to the corresponding operator from (2.1), X^ to 
Xb (see (1.24)) and vr^ ^ tt^ induces a Q^j— linear homomorphism from 

TlrC^ to the algebra of linear endomorphisms of Q^JX]. This - 
module is faithful and remains faithful when g, t take any complex 
values assuming that g 7^ is not a root of unity. It will be called the 
polynomial representation; the notation is 

v^'q;,[x,] = q;,[x,,6gs]. 

The images of the Yf, are called the difference-trigonometric Dunkl 
operators. 
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The polynomial representation is the —module induced from the 
one-dimensional representation Tj i— tj, Yj, q'^iPk^b) ^j^g affine 
Hecke subalgebra TYy = {Ti^Y^j). Here we extend the ring of constants 

to q;,,. 



2.1. Macdonald polynomials. There are two equivalent definitions 
of the nonsymmetric Macdonald polynomials ^ denoted by Ef,{X) = E^'^ 
for b & B; they belong to Q(g, t)'[Xa,a G B]. The first definition is 
based on the truncated theta function due to Macdonald: 

[l -X^qi)il-X-^qi-^') 



(2,3) MA-;t) =m"I(X) = n IfTT- V . _ v-u 



We will mainly consider Laurent series with the coefficients 

in the ring Q[ty] [[g^]] for u E ur = {usht, t'lng}- The constant term of a 
Laurent series f{X) will be denoted by (/). Then 

(2-4) (/i) = J_J_ _t^g(Pfe,a)+ji'a)(l _t-lg(pfc,a)+J>c.)- 

Recall that = g^"^'" = g^". This equality is equivalent to the 
Macdonald constant term conjecture proved in complete generality in 
[C2]. 

Let /io = /i/ (/i). The coefficients of the Laurent series /io are from 
the ring Q^^^. 

The polynomials Ei, are uniquely determined from the relations 
(2.5) ^6 - e ©e^bQ;,tX„ (EfeX- Vo) = for B3c^h. 

for generic g, t and form a basis in Qg^[Xb]. 

This definition is due to Macdonald (for /cght = ^ing G Z+), who 
extended the construction from [Op]. The general (reduced) case was 
considered in [C4]. 

Another approach is based on the F— operators. We continue using 
the same notation X, F, T for these operators acting in the polynomial 
representation. Let Xa{q'^) = g^"'''^ for a,b E P. 



Proposition 2.1. The polynomials {Eh,b G B} are unique (up to 
proportionality) eigenfunctions of the operators Ya (a G P) acting in 
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Q'jx] ■■ 

(2.6) YaiE,) = q(-+'P'')-i-'h) E, for 6„ = 6 - ul\pk), 
Uh = n^^h is from Proposition 1.1, = 7r6((— p^)). 

□ 

The coefficients of the Macdonald polynomials are well known to be 
rational functions in terms of qu,t^. In this paper, we use that these 
coefficients actually belong to Q^^^, i.e., are well defined when = 
for all It readily follows from (2.17) below. 

2.2. Symmetric polynomials. Following Proposition 2.1, the sym- 

(k) 

metric Macdonald polynomials Pb = P^ can be introduced as eigen- 
functions of the M^— invariant difference operators 

(2.7) V=RedH/( Yl ya') ioT a+eB+, 

a'£W{a+) 

where Redv^ is the restriction to the space V ^ of Vl^— invariants of V. 
Explicitly, 

La4Pb_) = q^''^^'''\ Yl q~^'''^'-~"'^)Pb-, b^EB^, 

a'eWia+) 

(2.8) n_= Yl 

mod ®c-yb- Q{q,t)Xc. 

beW{b^) 

These polynomials were introduced in [M2, Ml]. They were used for 
the ffist time in Kadell's unpublished work (classical root systems). In 
the case of Ai, they are due to Rogers. 

The connection between E and P is as follows 

Pb_ = Pb+Eb^, 6- e 5„, 6+ = wo{b.), 

(2.9) P;,^= Y ^^^^^ 

cgVF(f)+) 

Wc E W is the element of the least length such that c = Wc{b^). Taking 
the complete t— symmetrization P here (with the summation over all 
w), one obtains Pb_ up to proportionality. See [Op, M3, C4]. 

There are two different kinds of inner products in V from [C8] and 
other works. In the symmetric setting, they essentially coincide. We 
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will need here only the inner products of the symmetric polynomials 
Ph for b = b_ : 

(2.10) (A(X)P,(X-i)/io) 



a>0 j=0 



'^'''=11 n V (l-qiX^{qP.)){l-qi+^X^{qP>^)) 



2.3. Using intertwiners. The following map can be uniquely ex- 
tended to an automorphism of TH!' where proper fractional powers 
of q are added (see [C1],[C4],[C7]): 

r+ : Xii^ Xb, TTj. 1—^ q 2 X^Hr, 1— s> X^Y^q 2 , 

(2.11) r+ : To X^^{t^T^^}, Y^ ^ X^'{toT,-'}T,^. 

This automorphism fixes Tj {i > 1), t,y, q and fractional powers of q. 

The y— intertwiners serve as creation operators in the theory of non- 
symmetric Macdonald polynomials. Following [C6, C8], let 

(2.12) = T+{T,) + {U - l)(^a.(g'«) - 1)~\0 < z < n. 

We will use the pairing from (1.8) and the affine action w((c)) from 
(1.7). 

Theorem 2.2. Given c & B, < i < n such that (oj, c + d) > 0, 

(2.13) g(c,c)/2-(fe,6)/2^^ ^ ^c^^^) b=s4c)). 

If {ai,c + d) = 0, then 

(2.14) r+(ri)(Ej = 0<t<n, 

which results in the relations Si{Ec) = as i > 0. For b = 7ij-{{c)), 
where the indices r are from O', 

(2.15) qic,c)/2-m/2j^^ ^ r+(7r,)(i?e) = X^,.q-^'^-^^^/^7Tr{E,). 
Also T+{7ir){Ec) 7^ Ec for TTr j^id, sincc TTr{{cj) 7^ c for any c E B. 

□ 
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If (ai,c) > and i > 0, then the set X^iib) is obtained from A(-7rc) 
by adding [a,(c_,a)] for a = Uc(tti) G R- and (c_,a^) = {c,a^) > 
0. When i = and (oq, c + d) = — (c, + 1 > 0, then the root 
[a, (c_,q;) + 1] is added to A(7rc) for a = Uc{—'&) = G R- and 
(c_,a) = -(c,i9) > 0. 

In each of these two cases, (Q;j,M~^(p)) = (a,p) < and the powers 
of ty in 

(2.16) ^a,(g'«) = = 

V 

are from Z+ with that of ti strictly positive. 

Due to Theorem 2.2 (see also [C6], Corollary 5.3), the polynomial 
E\, exists if 

(2.17) n (l-?^^«(/'^))7^0. 

[-a, j]GA(7rb) 

If 6 G -B- and the latter inequality holds for 6+ = wq(\)) G 5+, then the 
symmetric polynomial Ph is well defined. If = for all i^, then E— 
polynomials and P— poljTiomials are always well defined, which gives 
that their coefficients are polynomials in terms of q. 

2 A. Spherical polynomials. The following renormalization of the i?- 
polynomials is of major importance in the Fourier analysis (see [C4]): 

(2.18) Sb = Eb{X){Eb{q-P''))-\ where 6 G fi. 



qLX^iqP") 



[a,j]e\' (TTt) 

This definition requires the t— localization. 

We call them nonsymmetric spherical polynomials . Formula (2.18) 
is the Macdonald evaluation conjecture in the nonsymmetric variant 
from [C4]. See [C3] for the symmetric evaluation conjecture. 

The following duality formula holds for b,c ^ B : 

(2.19) £b{q'») = £,iq'»), b^ = b-u^\p,), 

which is the main justification of the definition of Sf,. 

Given b & B, the polynomial Sb is well defined for g, t G C* if 

(2.20) n - QitMq'')) + 0. 

[aj]6A' (TTb) 
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In the symmetric setting, 

(2.21) Vt = P„(X)(Pi(q-'"'))-^ where 6 e B_ , 

The symmetric duahty reads as follows: 

(2.22) Vb{q'-'') = Vc{q'-''), for b,ceB.. 

The norm formula becomes entirely conceptual: 

(2.23) ((n(X)n(X-i)/io))-^ = Yl /^K)/i(id)-\ 

where fi{w) = f^iw(lq~^''))) for w G ly. 

It is a direct corollary of the fact that the Fourier transform sends the 
P— polynomials to the delta-functions; see [C8]. 

2.5. The limit t ^ 0. Let by the reduction of m!' for = 0, 
where u G z/_r. It can be called the nil-DAHA or the crystal DAHA. 
The polynomials E^, P^. are well defined and linearly generate V and 
correspondingly; V = Qq[Xi,,h G P], where Qg = Q(g^/™) is the 
ring of polynomials in terms of g^^™ with m from the definition of j. 
We will denote Ti{t = 0) by T^. 

Theorem 2.2 holds under this specialization and gives quite a con- 
structive approach to the P— polynomials. The intertwiners from 
(2.12) that appear in the formulas for Eb are all in the form r+(Ti) + 1 
in this limit. It is directly connected with the fact that T/ = Tj + 1 
satisfy the same homogeneous Coxeter relations as {Tj, < i < n} do, 
a special feature of the nil-DAHA. It readily results from the theory of 
intertwiners, but, of course, can be checked directly too. 

The action of vr^ on {T/} by conjugation obviously remains un- 
changed. Thus relations (i,ii) from Definition 1.2 hold and, given 
w G W, the element = T^rTl^ ■ ■ ■ T/^ does not depend on the choice 
of the reduced decomposition w = UrS^ ■ ■ - Si^. For instance, opera- 
tors = r+(Ti^J for i = l,...,n are pairwise commutative and, 
importantly, VT— invariant. 
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Indeed, one has: U'^ = Yli=ii^'i)"' -S- ^ b = —J^^i^i- P^O" 
vided that all rij > 0, the reduced decomposition b = = uio^6+ holds 
for the longest element Wq ^ W and 6+ = Wo{b) G 5+; see (1.18). Thus 
is divisible on the left by (Tj + 1) for any i > and therefore divis- 
ible by the VT— symmetrizer on the left. It results in the ly— invariance 
of Pfe for any b G B_. 

The ly— invariance of {life, b G B_} simplifies significantly the rela- 
tion of the i?— polynomials to the P— polynomials: 

(2.24) Pfo = Eb for b = b^eB^. 

In more detail, we have the following explicit proposition. 

Proposition 2.3. (i) In the representation V of TH^ , the polynomial 
r+(Ti)(l) equals q''^ Eb for w = 7rb,beB, Vb e Q. 
(a) In the symmetric case, 

(2.25) n;(l) = q'-'Pb for b e B^ , G Q, 

where II- can be replaced by their restrictions RedwO^'i) to , which 
are pairwise commutative W— invariant difference operators. □ 

It is important that only positive powers of q appear in the coeffi- 
cients of Eb- The coefficients of these g— polynomials are non-negative 
integers. One can obtain it from the interpretation via Demazure char- 
acters or using the intertwiners (we are going to discuss it in further 
papers). As g ^ 0, the polynomials Pb_ become the classical finite di- 
mensional Lie characters, which can be seen, for instance, from (2.29) 
below. 

For the affine root systems considered in this paper (with in terms 
of the maximal short root d) , the connection was established between 
the polynomials Eb{t oo) and the Demazure characters of the cor- 
responding irreducible affine Lie algebras. See [San] and, especially, 
[lonl]. Theorem 1. Paper [lonl] is based on the technique of intertwin- 
ers (from [KS] in the GL^— case and [C6] for arbitrary reduced root 
systems). We will not discuss this important relation in this paper. 

There is a relation between the limit t — > used here and the one 
t ^ oo. It goes through the general formula 

(2.26) El = n ^^^"''^"'''^""^^-ol^.w), where 

V&Vr 

X* = X~\ q* = q-\ t* = t-\ c^ib) = -woib), 



WHITTAKER LIMITS OF DIFFERENCE SPHERICAL FUNCTIONS 23 

form [C8] and other author's works. This connection is especially sim- 
ple for the symmetric polynomials: Pb{X)* = Pij{X~^) for b = 6_, 
i.e., Pb = Pb{t ^ 0) = P,(6)(t ^ oo). We use that Pb{X~^) = P,{b){X). 
This connects the g— Hermite polynomials and the Demazure characters 
for b = b^. 

Concerning the orthogonality of P, the denominator of the /i— function 
from (2.27) vanishes in the limit: 

oo 

(2.27) JI = n 11(1 - X^qi){l - X-'qi+^). 
The constant term formula becomes a well-known identity: 

n oo ^ 

(2.28) (/I) = -, where = q'^ 

i=l j=l ^ 1i 

For 6, c G -B- , the norm formula from (2.10) reads as: 

(2.29) {Pb{X)nX'')K) = 5b.ll n (1-?^)- 

i=i j=i 



3. Spherical and Whittaker functions 

We will begin with the identities involving the Gaussians, which are 
essentially from [C5]; then their limits t — > will be considered. 

The second part of this section is devoted to the Whittaker limit of 
the q, t— spherical function from [C5], which results in a formula for the 
g— Whittaker function in terms of the P— polynomials. 

We note that the Whittaker limit is a general procedure that can be 
applied to any solutions of the Macdonald eigenvalue problem (and its 
various degenerations and generalizations). 

3.1. Gauss-type integrals. By the Gaussians ^ we mean 

(3.1) r = Y. ^'^''""^^^ 7^ = E 

beB beB 

The multiplication by 7® preserves the space of Laurent series with 
coefficients in Q[t][[g2ft]], where rh{B,B) = Z is from the definition 
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of Q'qt- Accordingly, the coefficients must be taken from Q[t] [[i?"^^]] 
when the Gaussian 7® is taken. 

We will also use the real Gaussians defined as 

(3.2) = g^^'/^ where X^, = q'^x^ = {x, b),x^ = a;„. x^v. 

i 

Note that considering 7®=© as holomorphic functions (provided that 
|g| > 1 and, respectively, |g| < 1) the functions 7®/7 and 7® 7 are 
5— periodic in terms of x. 

The q—Mehta—Macdonald identity from [C5] 

(3.3) = n n( . ) 

provides the normalization constant for the g— Gauss integrals 

(3.4) (n(X)P,(X)7®^o) 

= g^^'^-(''+^''"=)Pc(g'-^^)n(g-''''-)(7®/io), 

where h,c E B_. Obviously, it implies the duality formula (2.22). 
Formula (3.4) can be naturally extended to the i?— polynomials (the 
proof even becomes simpler), but we do not need it in this paper. 

There are counterparts of (3.4) for 7® (treated as an analytic func- 
tion for |g| > 1), and for the Jackson summation taken instead of the 
constant term functional. See [C5, C8]. The considerations from this 
paper can be readily extended to these cases. 

Taking the limit. Let us tend t — > in (3.4). The definition of the 
P— polynomials implies that 

(3.5) limg-('='''^)Pc(g""^^) = for c G P_ , c+ = ^^(c). 

Note that — (c, p^) = (c+,pfc). For instance, it matches the evaluation 
formula in (2.21): limj_o g"^''''"^^c(g"^') = 1- 
We come to the following formulas (c G P_ ): 

n 00 

(3.7) (P,(X)P,(X)7®7I„) = g^^'^X,,(g'')(7^/Io). 
Here X^^{q^) = = = Xb^(g'=). 
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3.2. Global spherical function. One of the main advantages of the 
technique of Gaussians is a possibility to introduce the spherical func- 
tion as a reproducing kernel of the Fourier transform from V7~^, the 
polynomial representation multiplied by the Gaussian 7"^, to the TU}' — 
module V7. We will need only the symmetric case here. We assume 
that |g| < 1, which makes the considerations "naturally" compatible 
with the limit t — 0. In this setting, the construction below is directly 
related to the identities (3.4) (correspondingly, (3.7) in the limit). 

We note that if the whole polynomial representation is considered, 
then the corresponding anti-involutions of generally, require the 
t— localizations. Correspondingly, the definition of the Whittaker limit 
of the nonsymmetric counterpart of formula (3.9) below (see [C5]) be- 
comes more subtle. 

We will use the notation 7a and 7a for the Gaussians defined for 
another set of variables A completely analogous to X (7x,7x are old 
7,7). Thus, 7a = 7(5''^) and 7a = '^{q^) = q^^^'^. We will also use 



Theorem 3.1. Provided that \q\ < 1, the function ^! from the relation 



is a well-defined Laurent series. It is an analytic function for all X, A 
and for any choice of t^ assuming that all P— polynomials exist ( the 
conditions \ty\ < 1 are sufficient). 

The function ^°(X, A) defined via (3.9) is meromorphic for all X, A 
and analytic apart from the zeros of^^^^ . Replacing 7® 7® lx^l\^ 
in this definition, the corresponding function will he denoted simply by 
^(X, A); it becomes totally analytic but not a (single-valued) function 
in terms of Xf,, A^. 

Both functions, *P°(X, A) and%^{X,A), are X ^ A— symmetric, W— 
invariant with respect to X and A and satisfy the following extension 
of the eigenvalue problem from (2.8): 



(3.10) K^mX.A)) = q^^'^'^'K A:/)^(X,A). 



(3.8) 




l'^lt^\X,A)/^^{qP-) 



(3.9) 




a'&W(a+) 
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□ 

We note that Pb{X-^)Pb{A) = P,biX)P,h{A~^) in (3.9); recall that 
^(x) = —wq{x) and P^b{X) = Pf,{X~^). Applying q to the summation 
index b does not change the result. Thus: 

A) = r(A,^) = r('^w,'^(A))- 

The following can be used for an abstract (i.e., without an explicit 
formula) definition of the function ^°(X, A). It goes through the spher- 
ical polynomials {Vc = Pc/ Pc{(l~ '''')■, c G 5-} with a common coefficient 
of proportionality: 

(3.11) r(x,,-) ^ ^ n n(r^^™ ■ 

Here we substitute A = cj, A = g'^" in the left-hand side of (3.10) 
and divide it by the Gaussian 7®. This formula can be considered as 
a g, t— generalization of the Shintani-Casselman-Shalika formula from 
[Shi, CS]. Its limit as t — > will be discussed in the next section. 

3.3. Global Whittaker function. We are now in a position to define 
the global q— Whittaker function *P°(X, A) from the relation 

(3.12) 7.^^:(X,A) ^ }^ \f^^,,^ rr^-^^A). 

Here we always assume that t^, — > for all v. The function is 
defined for 7"^ instead of 7®: 

(3.13) ^x(A:,A) = linig(^'^^-)<p(g-^'=X,A). 

More explicitly, provided that |g| < 1 (we will not show the dependence 
of A here and where it cannot lead to misunderstanding): 

(3.14) $.(X) ^ 7x lim (7."^^) (g^"'"') = 

K— >00 \ / 

lim g^^g- '"""'2'"'''' <p(g^-Pfe) = lim g^^'"^) <p(g^-^^). 

fc^oo fc— >oo 

In this definition, A remains untouched, so the limit is a IV— invariant 
function with respect to A. As a matter of fact, the key fact we need 
is the existence of the limit 

(3.15) lim ^(g^-''\A;g,t) = §(X,A;g) 

fc— >oo 
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for A;g, t) from (3.9). Let us calculate the Whittaker \& in full 
detail. It is essentially a generating function for the P— polynomials ; 
see Proposition 2.3 and related formulas for the definition of these 
polynomials. 

Theorem 3.2. (i) Provided that \q\ < 1, the Whittaker function is 
given by the formula 

where the power series in the right-hand side is well defined coefficient- 
wise and converges everywhere; see (3.5) and (2.29). The formula for 
^2,' is with ''ix^'^x^ instead of 7®7® and with the same summation in 
the right-hand side. 

(a) The ratio of the functions ^°(X, A), ^^.(X, A) is B— periodic 
with respect to X and A. The dependence on A is governed by (3.10) 

for the limits L^^ of the operators as t ^ upon X A : 

(3.17) it^ mix, A)) = X-i$.(X, A), = 

In terms of X , these functions satisfy the q—Toda system of difference 
equations: 

(3.18) Z,^(^.(X,A)) = ( Yl A-/)^.(X,A), 

(3.19) L,, ^hra g'^"-'"^) (g^^'^'^'Hr"^ ^a, r,Jg-(^'^'=)) , 
Tb{F{X))=F{q'X), r,X, = g(^''^)X,r, for 6 G C". 

Here the difference operators La_,_(a+ G 5+) from (2.8) are conjugated 
by the translation F-p^, (it is {pk)' in the notation from (1.2)) and then 
by the operator of multiplication by q^^'p^), □ 

We note that Xf,^ Ph{A~^) in the summation for \& can be replaced by 
Xft~^Pfe(A). Recall that, generally, Pb{X-^) = P,(b){X) and b ^ <^{b) = 
—wo{b) does not change the coefficients in the summation from (3.16). 
The formulas for operators are simple to calculate for minuscule 
a+; see [Et] (An) and the rank one case below for examples. 
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Concerning the notation, one can introduce the functions 
using the Whittaker hmits with A, A instead of X, x, but we do not 
need these functions in the paper. Nevertheless, we put x in (not 
always) to emphasize that the Whittaker limit makes the dependence 
on X and A asymmetric. 

The construction of the Toda operators in terms of the Macdonald 
operators (and their various degenerations) is essentially due to In- 
ozemtsev and Etingof. The paper [Et] contains a systematic consid- 
eration of the difference case.. This paper is mainly about GLn, but 
our (3.19) is quite analogous to the limiting procedure there, as was 
expected in Remark 1 at the end of [Et]. 

We remark that our g— Toda operators are "dual" to those from [Et, 
GLOl] (the translation operators must be replaced by their inverses), 
which is connected with our choice of the limit t ^ versus t — > oo in 
these papers. The relation will be discussed below in greater detail. 

Theorem 3.3. Continuing the previous theorem, let X = q'^ for c G 
. Then the Shintani-type identity holds: 

n oo ^ 

(3.20) 7^(l)fc(g^A) = PciA)l[l[{-—-) 

i=l j=l ^ 1i 

w/iere 7®(1) = J^beB^^^^"^- More explicitly, 



n oo ^ 

(3.21) = 7^(A)Pc(A)nn(^ — l)- 

11 1 — a,- 



i=l j=l 

Proof. Due to (3.12), 

t-*o 7©(g~^fc) 

Applying the identity (3.11) for X transposed with A (the duality) 
inside 



(3.22) ^,B(g'm{q',A)= Ito le, (^° A)) ■ 



(3.23) ,^ifmf, A) ^ ( n n t4^) ■ 
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Recall that 

(3.24) (7)p, = (7)c-p. = T^lg^-^^-)?^"""^, 

where we use that c is from B; see (3.8). Hence, 

Moving q^'^'Pi'^ from (3.25) to the denominator and combining it with 
Pc{q~P''), we apply (3.5): 

limg-(^'''*)p,(g-'"=) = 1. 

Finally, we move q'^^^"^ from (3.25) to the left-hand side of (3.23) 
and observe that q'' /^7®(g'^) does not depend on c, so it equals 7®(1). 

□ 

We note that by making g = in (3.21), we arrive at the trivial 
identity ^^(A; 5 = 0) = Pfe(A; g = 0), where Pfe(A; g = 0) is the classical 
character for the dominant weight WQ{h). 

The p— adic limit g — (in this setting) transforms (3.11) to the clas- 
sical Shintani-Casselman-Shalika formulas. See [C8] concerning the p— 
adic degeneration of the DAHA theory (the limit g ^ oo is considered 
there). 



3.4. One-dimensional theory. We will begin with the explicit for- 
mula for the P— polynomials in the case of Ai. The formulas for the 
Rogers polynomials are well known as well as for their limits as t — > 0. 
Such limits are the q—Hermite polynomials introduced by Szego and 
considered in many works; see, e.g., [AST]. Let us re-establish the 
formulas we need for these polynomials using the (nonsymmetric) in- 
tertwining operators. 

Let a = ai = s = Si, u = uji = p; so a = 2uj and the standard 
invariant form is {nuo^muj) = nm/2. Similarly, 

X = X^ = g^ X{qn = g"/', r(F(X)) = F{q'/''X), 

i.e., x{nuj) = n/2, T{x) = x + 1/2, TX = q^/^XT. 

We will also use it = sT : X g^/^X~^; then vr^ =id and Y = = 
ttT in DAHA of type Ai. Concerning the Gaussians, ( 

2x2 and 7 = g(^'^)/2 = g^'; note that 7(g"'^) = g"'('^''^)/2 = g"V4_ Also, 
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X — pk = r~'^(x) = X — k/2 and g(^'^fc) = q^'^ in the formulas for the 
Whittaker hmit. Here t = = for k E C 
In the hmit t^O,T = T{t = 0) and 

P = P =E 

For instance, Pq = 1, Pi = X + X'^, 

P2 + ^ 1 + g, P3 = + X-3 + i^(x + X-^), 

1 - q 

(3.26) P4 = X^ + X-^ + ^(X^ + X-^) + ' ~ f ■ 

l-q (l-g)(l-g2) 

GeneraUy, for the monomial symmetric functions Mq = 1, M„ = X" + 
X"" for n > 1, 



[n/2] 



(3.27) Pn-M^^Y. - , "\ n 



The norm formulas from (2.27), (2.28), (2.29) read as follows: 



(3.28) (P^(X)P„(X)7I„) = 5^nX{{l-q'), 

i=i 

where m, n = 0, 1, . . ., 7*0 = V-l i'Po) the classical theta-function 
00 00 _ 

(3.29) /I = 11(1 - XV)(1 - X-'q^+'), (/I) = n 

j=0 j=l ^ 

Due to Theorem 2.2, the composition 7^ = (1 + T)Xtt is the raising 
operator for the P— polynomials. Namely, upon the restriction. Red, to 
the symmetric polynomials: 

ji^2p— 1 x~^r 

(3.30) q'^RiPn) = Pn+i, where R = Red(7^) = ^ _ . 

This readily gives (3.27). 

Rogers polynomials. The counterparts of these formulas for the 
Rogers polynomials are well-known (see, e.g., [Al] and [C8], Chapter 
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2). Let us list them for the sake of completeness. First, 

(/i) = Constant Term (/.) = H (i _ t2gj)(i _ g.) " 
Switching to = yu/ (yu), 

(3.32) (P„(X)P„(X)/io) = SmnJlj^ 



g^+i)(l-tV) 



as m, n = 0, 1, 2, . . . . The explicit formulas are as follows: 

(3.33) P. = M„ + y rf ^^—^ , M„_2,. 

The L— operators. We will begin with the formula for the g— Toda 
operator from (3.19): 

(3.34) L = lini T^^ L Tk g"*^" ) = (1 - X~^)T + 
for FfcX = t'^/^XPfc and the well-known operator 

1 - 

r H 

1-X2 1-X- 
diagonalizable in terms of Rogers' polynomials. We will also use: 

(3.36) = -f-'L-f = q'^^{XT + X^\T~' -T)) for 7 = g"'. 
For the straightforward specialization of L at t = 0, one has: 

L = limL = (1 -x2)-ir + (1 -x-2)-ir-\ 

(3.37) = 7-^X7 = -q^/\X - X-^)-\T - T-'). 

The latter operator is proportional to the so-called Askey -Wilson di- 
vided difference operator, which serves as the shift operator in the the- 
ory of Rogers' polynomials (with any t) and the basic hypergeometric 
function. See [Al] and also [C8], Chapter 2. 
Its defining property is the relation 

(3.38) L^(P„) = -g^/^(g" - g-")P„-i, n = 1, 2 . . . . 



(3.35) L = L^ = Red(r + tY-') = —T + -——^r~\ 
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Let us give a convenient reference concerning (3. 30), (3. 38): [OS], for- 
mulas (20-25). 

3.5. Whittaker function for A^. Provided that |g| < 1, we can now 

introduce the Whittaker function from the relation: 



The function A) is actually the generating function for q— 

Hermite polynomials. It is directly connected with the, so-called, qua- 
dratic g— exponential function; see [Sus] formulas (26), (27) and the ref- 
erences there. Its interpretation as a g— Whittaker function (upon the 
multiplication by the Gaussians) does not seem to have been noticed, 
although the difference equation for "^{X, A) was certainly known (for- 
mula (19) ibid.). The one-dimensional Shintani-type formulas, (3.43) 
below and especially its q, t— generalization, seem new. Some related 
formulas like (3.4) can be deduced from known identities (at level of 
e^E'e); the multidimensional theory is new. 

The power series \E'(X, A) converges everywhere. The A— dependence 
(see (3.17)) readily follows from (3.38): 

(3.40) l!^ (§(X, A)) = X §(X, A), = L^{X ^ A). 

In terms of X, the function satisfies the 'j— twisted g— Toda equa- 
tion, which reads as follows: 

(3.41) L^(^(X,A)) = (A + A-i)^(X,A). 

The q, t— case. The formula for the global spherical function 
reads as follows: 




(3.42) 



r(^,A)7.^7f^/7^(g') = ^(^,A) 




We come to a variant of the basic hypergeometric function. 
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Shintani-type formula. Let us consider Theorem 3.3 in the Pl- 
ease; we plug in X = g"'^/^ for n = 0, 1, . . . , . Then 

oo ^ 

(3.43) q-'/'^{q-/',A) = ^^{A)R,(A)l[(^-—-y 



i=i 



Recall that 7®(q''^) = S^gz Here the left-hand side and the 

right-hand side coincide as Laurent series or as analytic functions. 

This formula becomes a trivial identity for q = 0, i.e., in the case of 
the classical characters 



Pn{X;q = 0) 



X -X 



-1 



3.6. The case \q\ > 1. Generally, the Whittaker-type limiting pro- 
cedure as t ^ CX3 is naturally connected with the theory at |g| > 1 and 
can lead to new formulas. However, in the symmetric setting of this 
paper, there is a direct connection between the Whittaker functions 
defined for |g| < l,t ^ and |g| > l,t ^ oo, which we are going to 
discuss now. 

We follow [C5] and use 7® instead of 7® and 7 instead of 7"^. In 
the nonsymmetric setting, the corresponding global spherical function 
is really different from that for |g| < 1. However, there exists a simple 
connection in the symmetric case. 

The g, t— definition we need is as follows (cf. (3.9)): 



Pt(X) P„{\) 



(3.44, ^ ^)_,^----^ V>(^)n(-Y-);„> - 

where g > 1 and satisfies the claims of Theorem 3.1. The Whit- 
taker limiting procedure requires here taking t ^ 00 for ensuring the 
convergence. 

The formulas are: 

(3.45) 7.^7.^^: = 7." Hm ^:(g--X) 

(3.46) ^ §.(X, A-q) = J2 """^ 



2 



beB- Ili=illjii ' \^ ~ Qi ^] 
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Cf. (3.12) and (3.16). Here 

(3.47) limPfe(A;g,t) = lim P,{A-';q~\r') = P,{A~';q~'). 

Therefore \1/*(X, A; q) simply coincides with \[^(X~^, A; q~^) in the no- 
tation from (3.16). 

We conclude that satisfies the eigenvalue problem 

(3.48) L:^mx,A)) = i Yl K')^*iX,A), 

(3.49) LI ^\im g-('^-^'=)(g-(^''"=)(r-^ Vrpjg^^''"^)), 

where rfe(F(X)) = F{q^X). Compare with (3.19); the conjugation 
by qi^'Pk) there is replaced by the conjugation by q~^^'Pt=\ Thus the 
operators L*^{X,q~^) generalize those considered in [Et, GL02]. For 
instance, in the one-dimensional case in the notation from (3.34): 

(3.50) L* = lim (g"^'^ T7^ L q^"" ) 

t~*oo 

= (1 -x^)r + r-^ 



4. Harmonic analysis topics 

The real integration or Jackson integration is, generally, necessary 
when the Gaussian 7"^ in the constructions above is replaced by 7. 
A typical example is as follows. Let us consider the DAHA-Fourier 
transform in terms of the constant term functional (or using the imag- 
inary integration) in the space of Laurent polynomials multiplied by 
7~^. Then the inverse transform will involve the Jackson (or real) 
integration and the proper choice of the Gaussian is 7 instead of 7"^. 

Such "switch" of the Gaussians is necessary algebraically due to the 
properties of the involution of DAHA that governs the Fourier trans- 
form. Correspondingly, the contour of integration, real or imaginary, 
must ensure the convergence, i.e., its choice is of analytic nature. The 
direction is real for 7 and imaginary for 7"^. It is of course for |g| < 1; 
if |g| > 1 then it must be the other way round. Generally, especially, 
in the absence of the Gaussians (for instance, in the Harish- Chandra 
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theory), the directions, real or imaginary, are selected to match the 
growth estimates for the spherical function, used as kernels of the cor- 
responding transforms. 

We establish such estimates in the real direction. The theory appears 
surprisingly "precise", although the results of the paper are far from 
being complete. Only the first term of the asymptotic expansion is 
obtained. We note that in our setting, the global spherical function is 
periodic in the imaginary direction, so the imaginary growth estimates 
are irrelevant. We stick to the Jackson integration, which is actually 
very similar to the "classical" case of real integration; the estimates we 
obtain serve both theories. 



4.1. Growth estimates. It is possible to evaluate the growth of the 
global q, t— spherical function ^{X, A; q, t) from Theorem 3.1 in the real 
directions. Let < q < 1, = q^" (or, simply, t = q^) for E C 
provided the existence of all spherical symmetric polynomials {Vb_}, 
equivalently, provided that the polynomial representation is semisim- 
ple and the radical of the evaluation pairing vanishes (see [C9]). The 
assumption ^k^, > —1/hy for the Coxeter numbers hy = 1 + (p, {O'Y), 
where 9' = 9,{} for z/ = i^ing, I'sht, is sufficient (but not necessary). 

For X G C", let = u{x) where n(3?(x)) is a unique vector belong- 
ing to the closure C+ = Yl^=i ^+ ^« standard positive nonaffine 
Weyl chamber (t+ = ^"^-^M>o^i- 

Given a p— sequence of vectors x' = {x[, . . . ,Xp} C and a p— 
sequence of positive integers n = {ni, . . . , Up}, we use the dot-notation 
n ■ x' for Yl^=i ^i^'j- 

Theorem 4.1. (i) For arbitrary x, A G C", ky G C, we set: 

(4.1) A; *) -etJJl'-tlg,..) ^ «^ 

Given a real p— sequence x', let the components of n tend to +oo in an 
arbitrary way provided that (n ■ x')+ G Then the limit 

(4.2) lim <p2(x + n-x',A;g,A;) 

n— ►oo ' 

exists if^{X)^ G C+; moreover, it depends only on A and is nonzero for 
all such A. Here we choose x to ensure that ^9(g(^+"''' )++('^)+ "/''=) ^ 
for any n. 
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Under the same constraints, consider %^'^{x + n ■ x', A + n ■ A'; g, k) 
for a real p— sequence \' satisfying (n- A')^ G C+ . Then the limit exists 
too and is an absolute nonzero constant depending only on q,k. 

(a) In the case of the Whittaker function we remove k from the 
formulas and replace x+ by —x: 

(4.3) $;(,,A;,)Sl!|3W)<p.,,.,,A.,,. 

Provided that n ■ x' G — ^+ (it was not needed in the q,t—case), the 
claims from (i) hold true for 

(4.4) lim ^:(x + n-x',A + n- A';g). 

n— >oo ' 

Here ^^{q'^'^) is nonzero at [\ + n ■ A')_|_ — (a; + n ■ x') ; we continue 
to assume that (n ■ x')_|_ G C+ and, correspondingly, either 3?(A)_|_ G C+ 
for X'= or (n ■ A')+ G C+ . □ 

Let us comment on the proof. The justification of (i) involves the 
analysis of the corresponding difference equations for ^ in the limit of 
large x and/or large A, but we use the explicit formulas too. We note 
that the asymptotic difference equations provide the asymptotic limit 
(the factor in the definition of from (4.3)) only up to a periodic 
function. So we need to use that both, and are meromorphic. 

Using the asymptotic differential or difference equations for the anal- 
ysis of the limiting behavior of solutions of the corresponding equa- 
tions is standard. Given x, A, we restrict ourselves with the bi-lattice 
{x + P+, A + P+} and evaluate the values of there step-by-step using 
the corresponding difference equations. 

For instance, we treat the L— operator from (3.35) in the case of Ai 
as a recurrence for calculating the value at a; -|- (n -|- l)uji in terms of 
the values at x -|- nuji and x + {n — l)uJi, where the coefficients tend to 
constants as — > oo. The stabilization is of exponential type, which 
simplifies the necessary estimates. In the A— space, we use the x ^ \— 
duality and the A— counterpart of the L— operator. We note that in 
the rank one case, there is the classical theory by Birkhoff, which was 
developed recently in several works. 

The multi-dimensional case is not very different, as far as Theorem 
4.1 is concerned. Similar problems were considered in the theory of 
Quantum Knizhnik-Zamolodchikov equations. By the way, the equiv- 
alence of the QAKZ and the eigenvalue problem under consideration 
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(see [C8]), generally, can be used here. Our approach is direct. For ar- 
bitrary root systems, the formulas for the L— operators are not explicit, 
but we need only the stabilization estimates for their coefficients. This 
approach becomes significantly more transparent in the nonsymmetric 
theory, where we can use directly the intertwining operators instead of 
the difference relations. The nonsymmetric global functions have vari- 
ous symmetries including the transformation formulas under the action 
of the intertwining operators in both, the x— space and the A— space. 

Part (ii) is obtained as a limit of (i). Taking x' and A' real vec- 
tors is, actually, insignificant in the theorem. Since and are 
27rz log(g)P'^— periodic in the imaginary direction, it suffices to impose 
the conditions from (i,ii) for their real parts only. 

We note, that due to the claim that the limits do not depend on 
the particular way the integers {rij} approach the infinity, one can try 
to use the Shintani-type formulas. It is assuming that the uniqueness 
of is known in the corresponding analytic class of solutions of the 
spherical eigenvalue problem satisfying the Shintani-type formulas, i.e., 
among the solutions that "go through" the P— polynomials. To avoid 
misunderstanding, let us emphasize that we do not claim or use such 
uniqueness in this paper. An approach to its justification we have in 
mind employs the symmetry x ^ A or the passage to the nonsymmetric 
theory. 

4.2. Exact asymptotic formulas. Let us obtain them in the (most 
important) case A' = 0. As a matter of fact, we do not need (more 
general but less exact) Theorem 4.1, although the limiting difference 
equations are used in the proof of the main lemma. Recall that we 
deal with the function that is given by an explicit formula; the asymp- 
totic difference equations are natural and convenient here but can be 
replaced by straightforward analysis of based on the asymptotic 
theory of Macdonald polynomials. It is what we are going to do in this 
section. 

The key ingredient is the inverse of the positive half of the fi— 
function, a direct q, t— counterpart of the celebrated Harish- Chandra 
c— function [HC]: 



oo 



1 - tqXaql 

1 — XaQa 



(4.5) 



{x;q,t)= n n 



a6-R+ j=0 
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Theorem 4.2. (i) Provided the conditions of part (i) of Theorem 4-1 
for X' = 0, including 3?(A)_|_ G 

(4.6) lim ^°^{x + n-x.\X;q,k) = g{q,t)a{q^+;q,t) for 



n-^oo 



a>0 j- 



where (/i) zs i/ie constant term of fi from (2.4)- 

(a) Correspondingly, imposing n-x' G — C+ one? the other conditions 
in the Whittaker case, 

(4.8) lim rt(^ + n-x',A;g) = (/l) a(g"+; g, 0) 

n oo ^ 

/n contrast to this formula, assuming that (n ■ A')^ G t/ie X— limit 
does not depend on x: 

(4.9) \im^^{x,X + n-y-q) = (jl) = TT TT -. 

n— >(X) J. A J. A 1 „J 

i=l j=l ^ 'ii 

The limit remains the same if we substitute x ^ x + n- x.' in (4-9) for 
x' such that 3?((n ■ A')+ - n ■ x') G C+. 

Proof It suffices to calculate 

(4.10) lim <p2(c - pfe, A; g, /c), where c e B_ , 

c+— >oo ' 

and by c+ oo, we mean that (aj,c+) cxo for all z = I,-- - ,n. 
Recall that c+ = u'o(c), where c is always from i?_ in this calculation. 
Using the definition and formula (3.11), 

(4.11) mc-p,,X;q,k) ^ >^ > 



7® (g'*'+ ~'^)79 (^qpk ) 



("2) X 7^ n nf '° 



'^^^ ^ ae-R+ i=l 
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The special value Pc{q~^'') is given by (2.21); it is the exponent q^P'"^'^ 
times the product term, which will be combined (in the limit of large 
c+) with the product from (4.12). The result is exactly g{q,t), the 
constant term of the symmetrization of /i from [Ml, M2, C2]. 

We note that (yu) was obtained in this calculation without any refer- 
ence to its "true" meaning as the constant term of /x. It is interesting 
but not very much surprising; in [C8] the norm-formula for Macdo- 
nald polynomials (including the constant term formula) was actually 
deduced from the evaluation formula. Something similar occurs here. 

Since c ^ B (actually c G -B-), we can remove it from the theta- 
functions 7®(g'^~'''=) and 7®(q''*'+~'^), the multiplicators are the same as 
for the Gaussians and q~(^+~'^^^/'^. It gives: 

(4 13) 7^(g-^^)7^(g^) ^ (,,^_,^) 

The factor q^^^P^) will cancel the same term from Pf.{q~^'^) (in the de- 
nominator). The remaining part of (i) is taking the limit 

lim g-("'^+)Pe(g^), 

which is a subject of the following lemma. 
Lemma 4.3. Provided that |g| < 1 and '^{x^) G €+, 
lim g-(^'^+)p,(g^) = a(g^+;g,t), 

c+— >oo 

where the limit is pointwise. 

Proof. In the multiplicative notations, = g^"" M'^)-^) = X^(c) 

for w{^{x)) G C+, i.e., this monomial is from the leading symmetric 
monomial function of the Pc{X). Its coefficient is 1 by construction. 
One can assume here that w = 1 due to the VT— invariance of Pc- Then 
X~^Pc will be a power series in terms of X^. for z = 1, ■ ■ ■ , n. 

Calculating the corresponding difference equations (in the limit of 
large c+) is the most direct way to identify its expansion with <j{X). 
It suffices to uses the leading terms of the L— operators serving the 
symmetric Macdonald polynomials calculated in [C2], Proposition 3.4. 
Then we observe that (j{X) is a solution of this system of equations. 
It gives the required since both are power series in terms of Xq,- with 
the constant term 1. □ 
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The lemma gives (4.6). The Whittaker variants from (ii) are its 
straightforward hmits; the condition G — must be imposed in 
(4.8) and no such conditions are necessary in (4.9). Obtaining these 
two hmits via the Shintani-type formulas (3.20) seems possible as well, 
however, it requires knowing that these formulas determine uniquely 
in a proper class of functions, which we do not claim in this paper. 

□ 

Lemma 4.3 is known for the Askey- Wilson polynomials (see, e.g., 
[Is]). The Laurent expansion of the rank one /i— function is very ex- 
plicit, so it is straightforward. Paper [FZ] contains a comprehensive 
discussion of the Ai— case. In paper [Ru], the claim of the lemma was 
obtained in the A— case for the L^— convergence. It was conjectured 
there (with some explicit estimate) that the convergence is pointwise 
as well; see a discussion after formula (1.23). Paper [vDl] is an ex- 
tension of [Ru] to the case of arbitrary reduced root systems (for the 
strong L^— convergence). See also [vD2] for the case of the Koornwinder 
polynomials (the root system C^C„). 

Our operator approach (based on the asymptotic difference opera- 
tors) gives the pointwise convergence. We can, generally, answer Rui- 
jsenaars' question concerning the pointwise estimates in compact sets. 
However, we will not touch upon this (important) direction in this 
paper. 

As for Theorem 4.2, we think that its one-dimensional versions (for 
the basic hypergeometric function or its variants) are likely to be known. 



4.3. The Harish-Chandra formula. The corollary is an exact gen- 
eralization of the Harish-Chandra fundamental asymptotic formula for 
the classical spherical functions. Indeed, for x approaching oo in the 
directions x' (admissible in the sense of Theorem 4.1), asymptotically. 



(4. 14) lim (x + n ■ x'. A; q, k) 



n^oo 



Up to a simple py— invariant and 5— periodic factor C(x, A), depending 
of course on g, k (it is Z— periodic in terms of fc), we can switch to ^ 
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here, replacing all 7®(q'^) by 7 ^(g^) = q ^'^/'^. It gives that in the limit 
of large 3fJ(x)+ G C+, 

(4.15) A; g, A;) ~ A) ^(g, t) A+-Pfe)+(A+,p,) 

Corollary 4.4. VFe continue to assume that all spherical polynomial 
{^b-} exist; for instance, the conditions ^ —1/hy — Q+ for the 
Coxeter numbers hy of R are sufficient. Provided that 9ft(A)+ G 
the global spherical function A;g, /c) is bounded in terms of x as 

3 ^ 00 if and only if 

< (3?(A)+, a^) < ^{ki) for i = 1, . . . ,n, which implies ^{k^) > 0. 

//3?(A)+ G C+ allowed, then ^^(x, A; g, fc) asymptotically approaches 
a polynomial in terms of {xi} of degree no greater than n, the rank of 
the root system. □ 

The dependence of x in the right-hand side of (4.15) is as in the 
Harish- Chandra formula [HC]. Accordingly, Corollary 4.4 is a q,t— 
version of the description of the bounded spherical functions from [HJ]. 

The corresponding degeneration of ?ti (and all related objects) is 
the procedure g — *■ 1, where we set Xf, = e~^'' and Zh, Xb, k are con- 
sidered the basic new variables upon the degeneration. We take —Zb 
here because the base g is smaller than 1. The limit of the right-hand 
side of (4.15) can be readily controlled using the functional equation 
for the theta- function 7® (see below). Up to some renormalization, it 
becomes (for large 3fJ(z_|_)): 

Con.t n - n Flirt) Z = ^ = A = A,. 



i=l aeR^ 

The factor q^(^+'Pk\ which ensures the X ^ A— duality of the q,t— 
formula, vanishes in the limit. The duality collapses under the degen- 
eration to the Harish-Chandra theory, however, the evaluation formula 
survives. 

Technically, (4.15) matches the growth estimates for complex Lie 
groups because real Lie group result in the terms like r(A^/2) in this 
formula, which is not the case. 

The g — ^ 1 limit of the global spherical function is convenient to 
describe in a somewhat different normalization. Using the notations 
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from Corollary 4.2, we set 

def ^°{x,X-q,k) 7®(g^+)7®(l) 



(4.16) mx,X;q,k) 



g{q, t)a{q^+;q, t) ^e{q^+-Pk)^e{qPk] 



Cf. (4.1). Then '^'^{x,b — pk;q,k) = Pb{q^;q,k) for all symmetric 
Macdonald polynomials Pb, b G B_ . 

Apart from the zeros of 7®((3'^), this function is well defined for any 
A if all Macdonald polynomials {Pb,b G B^} exist. This condition is 
weaker than the existence of all spherical polynomials {Vb,b G B^} 
we imposed above; see (2.17)and (2.20). Moreover, if q'^ is not in the 
form for w G W, then is well defined for arbitrary k (i.e., 

the conditions for k necessary for the existence of {Pb, b G B_} are not 
needed) . 

Theorem 4.5. Provided that the Jack-Heckman-Opdam polynomials 
P'i^{z,X;k) = lim^^i P;,(e~^, g'^; g, g^) are well defined for all b G 5_ 
(a condition for k), given arbitrary complex z,X, the following limit 
exists : 

(4.17) ^'{z,X-k) = hm rt(e-^g^g,g'). 

This function is a W— invariant solution of the system of differential 
equations from [HO] and satisfies the following conditions : 

^',{z,-b-pk;k) = P'b{z,X;k) for beB^. 

Moreover, if X ^ W{B^ + p^), then the limit exists for any k. □ 

Here one can take complex q = exp(— 1/a + «0) provided that a > 
0, a — >■ cxD and —C/a < (p < C/a for a certain constant C. Numerical 
experiments show that here C can be arbitrarily large for any given 
(admissible) z, X, k, but we cannot justify it. 

The growth estimates for read as follows: 

(4.18) A; ~ lefa.,-..) yV.) ■ 

where the asymptotic equivalence must be understood as in Theorem 
4.1 under the conditions from (i) imposed there. The estimates are the 
simplest for such normalization, since "goes through" the Macdo- 
nald polynomials. 
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We note that "extending" (4.16) and (4.17) from the points in the 
form — 6_ — pk to all x is, generally, a non-trivial problem. We involve 
the growth estimates. 

The estimate (4.16) becomes exactly q~'y^+<^+~pk) up to a periodic 
function. The latter can be readily evaluated using the following (clas- 
sical) functional equation, a progenitor of the quantum Langlands cor- 
respondence. The following is a variant of the formulas that can be 
found in [Kac]. 

Lemma 4.6. Let A he the lattice dual to B with respect to the standard 
pairing {, ) in M", [B : A] the index from the theory of lattices; for 
instance, A = if B = Q and [B : A] = |P^/(5|~^. Then, picking 
u E C such that < < oo, 

(4.19) f/"'/2 ^ Xfef/'''/' = {V2^T^/[BTT] ^YaV'/^ 

beB aeA 

setting: U = exp{--), V = exp{--) forv = -\-,y = -^, 

where Xb = Ya = V^"- ; for complex u,v, we set = exp{—z/u) 
and = exp{—z/v). 

Recall that {x, x) /2 = x\ — xiX2 + X2 for Xi = x^. as B = P in the 
case of A2; correspondingly, {y,y)/2 = {yf + yiy2 + yl)/?) for A = Q, 
Hi = Uai = {2xi — Xi>)/ {2mu) , where i' = 3 — i, i = 1,2 for A^. □ 

Claim (i) of Theorem 4.1 can be naturally modified toward the Whit- 
taker limiting procedure as follows. 



4.4. When h — > 00 . Let us reformulate (4.21) entirely in terms of the 
function \E' from (3.9). Namely, provided the conditions from Theorem 
4.7,(i), the limit of the function 

(4.20) ^^[x,\-q,k) ^ {^^{q'^+'^-^>'))-'^iX,A;q,t) 

exists. Similarly, ^^{x,X;q) ={^^{q^+-'')y^ ^{X,A;q). The Whit- 
taker limit becomes simply: 

hm V[/(g-P^g^;g,g^-) = $(g^g^g). 

fc— »oo 

See (3.15). 



44 



IVAN CHEREDNIK 



Given real k'^^ > 0, let us replace k by k + n'k' for n' G N in (4.2) and 
analyze the limit 

(4.21) lim ^^{x + n-x',X + n- X';q,k + nk'). 

{n,n'}— >oo 

In the non-simply-laced case, n' can be treated as a 2— vector {n'^^^, n[^g} 
and n'-k' considered instead of n'k'; then both components are supposed 
to approach infinity (in this paper). 

Theorem 4.7. We represent k'^ = Uy + for non-negative real u^, 
and pick the directions x', A' such that 

(4.22) (a) (n-x)+-nV„ G C+ 3 3?(A)+ when A' = or 

(6) (n • x')+ — n' pu G C+ 9 (n ■ A')+ — n' when A' 7^ 

for all n, n' . Then the limit (4-21) exists subject to conditions from part 
(i) of Theorem 4-i, including the strict positivity requirement 3?(A)+ G 
£+. It does not depend on x in case (a) and is a x, X— constant under 
(b). If k' > then the limit does not depend on k too, i.e., depends 
only on x for (a) and is an absolute constant for (b). □ 

The justifications are based on the formulas for the asymptotic dif- 
ference equations for the functions under consideration. Theorem 4.1 
corresponds to the case k' = 0; then the limit does depend on k. The 
rule here is that the limit does not depend on the vectors x, X ot k 
involved in the limit, provided that the corresponding directions and 
the values of the vectors which are fixed are generic. 

The Whittaker limiting procedure can be treated as an extreme case 
of the theorem as follows. Let k = n'k' assuming that k' > and 
n = {n'}. We take A' = 0, x' = —pk'- Then the limit (4.21) still exists 
but now it depends on x (and depends on A too because we set A' = 0). 
Explicitly, 

^t(a;,A;g,A;) = (7Q(g^+-^))"^ ^(g^+"'^', g, g"''^'), 

since (x + n'x')+ = pk — x for sufficiently large n' . Actually, we do not 
need ^'■j- here; the correction factor (7®(g'^+~^))~^ does not depend on 
n' . We arrive at the procedure from (3.12). 

We believe that the following calculation is clarifying. Let us take 
generic extreme x' and A' in (4.22): 

x' = Pu, X' = py, so + A'+ - = 0. 
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Similar to the Whittaker case, we do not need ^'-j- here. Assuming that 
all and v^, are nonzero, 

lim ^(x + n' x', \ + n' A'; q,k + n'k') 



n'—ioo 



(4.23) = 



2 



beB+ Yli=l Yl\=i ~ 1i. 



for A'^ = wo{\~^). Thus, we obtain a non-constant dependence on x 
and A here, but the output is (one of the variants of) the multi-variable 
g— exponential function, i.e., significantly simpler than the Whittaker 
function. 

A Whittaker variant of this calculation is actually an extreme case 
of formula (4.9). It is: 

(4.24) lim ^(x (n ■ y')+, A + n ■ y'; g) 



(b,b) q 



xi,—X 



2 



beB+ lli=lllj=l 



n:Lin&/^(i-?r 



where we use the same y' for x and A (but in somewhat different 
way), assuming that 3?(n ■ y')+ G C+. Note the sign of (n ■ y')+; 
the growth estimates for the g— Whittaker functions considered above 
required taking the direction from the negative Weyl chamber. The 
proof is simple; we only need to know the leading coefficient of Pb is 1. 



Discussion. The theorems guarantee exponential growth (to be 
exact, no greater) of the function including the boundaries of the 
domains in the theorems. 

In more detail, the Gaussian-type corrections used in the definitions 
of ^—functions and the corresponding functions are not sufficient 
to ensure the existence of the limits on the boundary of the domains 
considered in Theorem 4.1 and 4.7. Even if they are sufficient for the 
convergence (as in the Whittaker case), then the limits can depend on 
the initial x, A. For instance, when (n ■ x') + , (n ■ A')^ belong to faces 
of the Weyl chamber the limits are expected to be connected with 
the spherical (and Whittaker) functions for subsystems of R. 

The role of the condition 3?(A)+ G as A' = is also important 
and not clarified in full. As it was claimed, if 3?(A)+ ^ then (4.20), 
generally, diverges, but the growth is polynomial. 



46 



IVAN CHEREDNIK 



A description of such and similar extreme situations and the corre- 
sponding asymptotic systems of difference equations is a natural chal- 
lenge. 

Numerical experiments in the rank one case confirm that the con- 
vergence condition (4.22) is sharp. It is not clear what happens if 
x' is taken non-proportional to pk' (especially in the non-simply-laced 
case when k = {/csht? ^ing})- Generally, for any x' G £+, the conver- 
gence of \l/(x + n'x', \;q,k + n'k')-^ is granted for < /c' < /Cq, where 
ko = ko{x') > 0. What is the formula for ko{x') and for which x' the 
limit exists at such extreme ko{x') ? 



4.5. Jackson integrals. We are going to integrate the product of 
two global spherical functions for the //—measure twisted by the plus- 
Gaussian. The previous section guarantees that the growth of this 
function in real directions is no greater than exponential. Due to the 
presence of the Gaussian, this is sufficient to ensure the convergence 
of the Jackson summations in the theorem below. This theorem is not 
from [C5], but its proof is based on the same technique (see also [C8]). 
Let us fix ^ G C" and define the Jackson summation as follows: 

{f)i = \W\-' where «;(0 + 6 = (M((0)- 

w£W,bGB 

Here the affine action of W from (1.7) is used; / can be any function 
well defined at the set {q'^^^^~^^}. Recall that the notation (/) was used 
for the constant term of a Laurent series /. We continue to assume 
that |g| < 1. 

As above, X„(g«) = q^°''^\'j{q'') = g(^'^)/2, [z,z) = EILi^i^"^' say, 
{z, z)/2 = zf — Z1Z2 + Z2 for A2. For instance, 

a6-B aG-B 

We will constantly use that (7)^ is periodic with respect to the sub- 
stitutions + b,b E B. As in [C5], let us introduce the function 

Ji{X;t) = fi~^{X;t~^), a counterpart of fi in the theory of Jackson 
integration; all ta must be replaced by t~^. 
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Let US set = Explicitly, usine the sets 

\{bw) =R+n {bw)-\-R+), 

„ //"1/2 _ .1/2 

(4.25) /..(g-(^)-^^) = n V. -V.' . = • 

[a,Vci]GA(6i«) y 
def P/ 



In terms of the action = f{q^ ^^'^'') on functions f{q^ 

(4.26) (//i.)^ 



mm) 

w&W,b&B 

Theorem 4.8. For arbitrary weights A = q'^, A' = q^' , 
(4.27) (7®(g^^))' (r A)r (^-\ A ') 7/i.)€ 

= hA..>£7?7§r(A.A') n iKri^rWi). 

ae/?+ i=i ''a 



|l^|/^(g«) n Hv i_ -(P.,"-)+i 



□ 



In these formulas, ti, are arbitrary provided the existence of all {Vb}- 
The products are considered as the limits if ki, G Z+ \ {0}. The nor- 
malization factor is obtained by taking A = A' = q~P''. Indeed, 

r(x,g--) = n n( ._ ':>.v,,, )r(/Ng---) and 

(r(^,g-''^)r(^,g"''^)7/i.)? 

(4.28) = n n(r^5TS^)rfe-,.-«){7^.)^ 

due to formula (3.11). 

Theorem 4.8 generalizes that from [C5] (the case of the Macdonald 
polynomials). The best way of obtaining the identities from (4.27) is via 
the interpretation of ^°(X, A) as the reproducing kernel of the Fourier 
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transform, but here the nonsymmetric setting is more convenient. We 
are going to follow this approach in the next paper (s). 

4.6. The special case ^ = —pk- The theory of Jackson-Gauss in- 
tegrals is essentially algebraic, similar to that for the constant term 
functional. Analytically, we need only the exponential growth of ^ 
in real directions; (4.2) is more than sufficient. The growth estimates 
can be equally used in the theory based on the real integration instead 
of the Jackson summation. This theory is a g— generalization of the 
so-called non-compact case in the harmonic analysis on the symmetric 
spaces. Formulas like (4.27) hold in such theory but the corresponding 
factors of proportionality (generally, periodic functions in terms of X 
and A) are not calculated so far with a reservation about the Ai— case 
(see [C8], Etingof's theorem). 

There is a special case when (4.27) becomes a straightforward alge- 
braic exercise; it occurs for ^ = —pk taken as the starting point of the 
Jackson summation. In this case, jj,{q'^^^^~^'') is nonzero if and only if 
bw = Tib = bu^^, i.e., at fej = TTb{{—Pk)) = b — u^^{pk) in the notations 
from Proposition 1.1. One has: 

(4.29) ,.(,'.) =.^"-"'n*t''"" n C'^dCZ^ ' 

where X^Hb) = { [a,j] \ [-a,Uaj] E \{7ib) }. Then (7)5 = (7)^^ and 

(4.30) (7^.)-„ = |Mr'(7>„- n Uir^pru^)- 

a£R+ j=l ^ '^a 

Formula (4.28) reads as follows: 

\w\ ii^iq^r (r(^,A)r(^^\A')7/i.)-p. 

(4.31) = (7)..7?7? n n( . "'r,...-,., ) r(A,A% 

a€R+ j=l ^ ''Q 

It is important to note that (4.31) is not a new identity. To be more 
precise, it formally results from the definition of the duality of the 
P— polynomials and the Shintani-type relations from (3.11). In a sense. 
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the Jackson integrals trivializes at the special ^ = — pfc, which is anal- 
ogous to the normalization condition in the theory of spherical func- 
tions. Actually relations (3.11) were deduced in [C5] from the general 
^—theory of Jackson integration, so this analogy is with reservations. 

4.7. Taking the limit. Let us interpret the identity (4.31) upon the 
Whittaker limit. The Jackson summation will be now over B: (/)<> = 
l^feGS fil'^)'^ notice that there is no factor versus the previous def- 
inition. For instance, (7)0 = 7®(1)- The corresponding /I— measure is 
nonzero only on B^: 

(4.32) Mi'^) = n n ^-iir'- 

i=i j=i 

We come to the following "reformulation" of the definition of 

(7)o(g(^'')r(^-\A')77io)o 

n 00 _ 

(4.33) = 7a^7.^ nn(T^) 

1=1 j=i ^ % 

Here the Shintani-type formulas were employed. 

It is instructional to obtain (4.33) as a Whittaker-type limit of (4.31). 
We suggest the following way. 

First, let us make k a positive integer, to be exact, k = N for = 
{Ni, G N}. Then G P+ and, for instance, {'j)x+p^ = {'y)\, which will 
be used constantly. Second, let us renormalize the /i— measure (4.29): 

/I.(g''«) =g-'(''-''"=V.(g'«) 

(4.34) =n*t'" n (f^^frWj)- 

The limit of Jl,{q^^) as — 00 exists for any b & B and is nonzero 
only for b = b^. Namely, 

hm M(l''-)=M<l'^)- 

Third, we will use the following property of the spherical polynomials: 

(4.35) lim n(g^"^^) = for b = 6_, 

A''—* 00 

which is a reformulation of (3.5). 
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Forth, we observe that the condition pn & B guarantees that 



b 



for any x. Therefore the hmiting procedure for obtaining from 
from (3.12) coincides with that for *P from (3.13): 

(4.36) ^°(X,A) = hm g(^'''iv)«p°(g-PivX, A). 



Replacing now A by A — pjy in (4.31), one obtains: 

\'\^,,, \ 

PN 



(4.37) (7)0 ' {^^{qP-)f{^°{X, q^-^-)^\X-\ A') 7/i.) 



2 00 , (pjv,o^)+j 

Qe-R+ j=l ''a y° 

Here is not present due to our definition of the Jackson summation 
in the Whittaker case. 

We can restrict ourselves only with 6 = since the other h appear 
in (4.34) with strictly positive t— factors Hi. ^j'"^"''''- Then the left-hand 
side of (4.37) modulo higher powers of t is as follows: LHS mod (t) = 



v2 



Q g^/I.(6_)<p°(g^-"''iv,g^-p^){g(''-'Piv)<p°(g-fe+-Piv,A')} 



for g = {g''-/'7®(g^-)/(7)p,}7®(/-) = 7®(/-) 

and n n 



aeR+ j=l ^ ''a 

Transforming correspondingly the right-hand side of (4.37), one ar- 
rives at 

RHS = q-^ {g"^(7)A}7A®nM?^''''-^r(9'"'-,A')}. 

The term q~pN/'^ can be moved to the LHS and combined with Q, 
namely, 

/-/'7^(g^-) = (7)p,. 
One n can be reduced in the LHS mod it) and the RHS. 
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Then we use (4.35) for 'P{q^ ^^'^ and the definition of the Whit- 
taker hmit (4.36) for ^°{q~^+~pN ^ X') and for q3°(g^-^iv , A'). Replacing 

(back) q {'^)x by 7® (A) and changing the summation set in the LHS 
from B_ to B^, we eventually obtain (4.33). 

This calculation is expected to be a sample for the general ^—Jackson 
integration theory in the Whittaker case (presumably, for the real inte- 
gration too); it will be discussed elsewhere. We note that the term 
in the integrand of (4.33) can be naturally combined with ^ = 
and "eliminated" upon the change of variables x + \ ^ x. However 
this substitution will change the summation set from 5+ to A + 
i.e., the general Jackson summation (with an arbitrary starting vector) 
naturally emerges even in the special case under consideration. 

The extreme case. There is no "natural" way to eliminate A, A' 
from (4.33) by evaluating this formula at certain special points. Gen- 
erally, such elimination is a standard way of discovering new identities 
that contain only q. Another possibility is in taking A ~ A' — oo for 
A = g'^. A' = q^\ t Let us perform this calculation in detail. 

We will use (4.9): 

(4.38) lim $(A + (n-y )+,A' + n-y ;g) 



(b,b) q^b ^'b 



b&B+ Y\i=lY]^j=l ' 0- ~ 



where (n ■ y')+ G C+. In this limit, formula (4.33) reads as: 



(4.39) = 



where we canceled out {(i) = YYi=i YYiLii^~^i) both sides. Moving 



beB+ Ili=l ltj=l \^ ~ ^i. 

b from the arguments of 7® and using that (7)0 = 7®(1), we come to 
an identical equality. No new formulas appear in this way. 

Discussion. We think that the growth estimates and formula (4.33) 
show great potential of the g— theory of Whittaker functions in har- 
monic analysis. For instance, an immediate interpretation of (4.33) is 
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the fact that the global g— Whittaker function multiplied by the Gauss- 
ian is essentially invariant with respect to the g— Fourier- Jackson trans- 
form for the measure JI^ from (4.32), which is very much standard in 
the theory of g— functions. 

This paper seems a convincing demonstration of the key role of 
Shintani-type formulas in the theory of spherical and Whittaker func- 
tions. Interestingly, quite a few analytic facts (e.g, the g— generahzation 
of the Harish-Chandra asymptotic formula) are directly related to these 
formulas. This is different from the differential setting and makes the 
g— theory significantly more algebraic than the classical harmonic anal- 
ysis on the symmetric spaces. 

We would like to mention that global spherical and Whittaker func- 
tions are expected to have properties similar to celebrated Ramanujan's 
mock theta functions, including the theory at |g| = 1 and certain (but 
not direct) counterparts of Maas forms. To be more exact, the natu- 
ral objects associated with g— spherical functions are Maas-type theta 
functions , which are not holomorphic in terms of x, A but satisfy the 
modular equation with respect to g. 

It must not be very surprising because the basic hypergeometric func- 
tion is known to be related to (some) mock functions. Our global sphe- 
rical functions are its multi-variable generalizations. 
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